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Abstract: Understanding the deflection of liglit by a massive deflector, 
as well as the associated gravitational lens phenomena, require the use of 
the theory of General Relativity. I consider here a classical analogy, based 
on Newton's equation of motion for massive particles. These particles are 
emitted by a distant source and deflected by the gravitational fleld of a 
(opaque) star or of a (transparent) galaxy. The dependence of the deviation 
angle D on the impact parameter 6, and the - Euclidean - geometry of the 
(source, deflector, earth) triplet, imply that different particle trajectories 
may reach an earth based observer. Since D{h) does not depend on the mass 
of the particles, a (Newtonian) flavor of gravitational lenses phenomena is 
naively obtained by setting the particles' velocity equal to the speed of light. 
Orders of magnitude are obtained through this classical approach, and are 
compared to the General Relativity results. 
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1 Introduction 



Gravitational lenses phenomena [H 12 El HI 13 E] rest on the gravitational 
deflection of light, and their explanation require the use of General Relativity 
(the gravitational deflection of light by the sun, of order 10~^ rd, was indeed 
a major issue at the very beginning of the theory). This theory uses a rather 
heavy technical machinery, which is not easily accessible to undergraduate 
students. One of the motivations of the present paper is to offer, through 
a classical approach, a feeling (and orders of magnitude) for gravitational 
lenses phenomena. 

Apart from extreme cases (neutron stars, black holes,...), the gravitational 
defection of light in General Relativity (GR) is small, and can be described 
by a weak gravity approximation jTll21E|- For quasistationary isolated mass 
distributions, the evolution of the unit tangent vector (e(s)) to a light ray 
(r(s)) is given to lowest order (see below) by 

where c is the speed of light, U the gravitational potential created by static 
mass distributions, and (y±U = VU — e^e-VU)). The weak gravity approx- 
imation corresponds to ^ << 1, and "lowest order" means that the typical 
velocity v of the mass distribution is small compared to c: A more rigorous 
calculation indeed shows that there appear other terms in the r.h.s of equa- 
tion (^, among which the lowest order term is of order 0(-)) (see eq. 4.17, 
p. 124 of 0). 

Equation (Q) bears some resemblance to a classical equation of motion; 
it is indeed a familiar remark that Newton's classical equation for a massive 
particle in the gravitational field of a deflector leads to a mass independent 
deflection angle (identity of inertial and gravitational masses). In this paper, 
we study the classical mechanical problem defined by the triplet (particle 
source (S), deflector (A), earth based observer (E)). An important ingredient 
of this study is the dependence of the gravitational deflection angle D on the 
impact parameter b of the particles. Since the (Euclidean) distances between 
(S), (A) and (E) are finite, the exact calculation of D{b) is rather involved, 
even if one takes advantage of the central character of the gravitational force 
and of some invariant properties (we use here the Runge-Lenz vector and 
tensor). A rather quantitative approach is used here to derive the main 
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features of the deflection D{b) for the case of spherical deflectors (mass M, 
radius R). 

To make contact with the (GR) weak gravity approximation, we consider 
a similar approximation for the classical mechanical study, namely we con- 
sider the limit (Gravitational energy << Kinetic energy). Further, if (foo) 
is a measure of the particle velocity as it leaves the source, one has at some 
point to make the (delicate) correspondence v^o ^ cin the classical mechan- 
ical problem. In this paper we follow this view of the gravitational deflection 
of light to obtain orders of magnitude, using typical undergraduate skills. It 
is shown that, in some cases, several particle trajectories may reach observer 
(E). These trajectories translate into a basic gravitational lensing effect: the 
source (S) may give multiple (photonic) signals for observer (E). The time 
delay between the reception of these different trajectories (or signals) is also 
considered. 

The plan of the paper is as follows. The gravitational action of the deflec- 
tor, such as the sun (opaque star) or the Milky Way (transparent galaxy), 
is studied in Section |21 For simplicity, we assume there that the source (S) 
and the observer (E) are both at infinity. When (S) and (E) are at a finite 
distance from the deflector, a little bit of geometry shows that that (S) may 
have multiple "images" (Sectional). Finally, an estimate of the time lag be- 
tween the reception of these "images" is obtained, via the analogy of particle 
mechanics and geometrical optics (Section 

2 Study of the classical mechanical deflection 
problem 

We consider a source (S) which emits, in an isotropic fashion, isovelocity 
particles (I^TI = Voo) of mass m. Unless otherwise stated, we will consider 
a point source. We focus our interest on the deflector and take the source 
and the observer at infinity. Let Voo = VooUi be the initial velocity of a 
test particle. The deflector (A) may be either opaque or transparent to the 
particles. The former case will be illustrated by the sun (radius Rq ~ 7-lO^m 
and mass Mq ~ 2 ■ lO^'^kg), and the latter by a spherical version of the Milky 
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Way ^ (radius Rq — 2 ■ 10^ light years (ly) ^ and mass Mq ^ 10^"^ Mq). 



2.1 The opaque detector 

The impact parameter b and deviation angle D are defined in Figure 1. The 
study of the classical Kepler problem for b > can be found in standard 
textbooks (see e.g. [7j). If (-Uj) and (uf) are the unit vectors along initial 
and final velocities, we have cos D{b) = Ui ■ Uf. 

The Runge-Lenz vector of the Kepler problem reads 

A = vxL-gMm- (2) 

r 

where r is the particule position measured from the center of the defiector 
(A), L = fxmv its angular momentum, and Q the universal gravity constant. 
As is easily checked, A is an invariant vector of the motion and writing 
Ai = Af ^ leads to 

tan — = T-TT (3) 
2 bvl ^ ^ 

The weak gravity condition means that the potential energy mil (b) is a small 
part of the total energy, that is << 1. The deviation is then given by 



The (GR) result is obtained from equation ([T)) as DcRib) ~ ^£3^. Apply- 
ing these results to the sun at grazing incidence yields Dgr{Rq) ~ 8 10~^ 
rd. 

Two remarks are in order (i) The naive identification {voo — >■ c), yields 
an expression which is half the exact (GR) result (ii) In classical mechanics, 
energy conservation implies that v'^{r) = v'^ + If one identifies foo with 

^see e.g. 



|http: / /hypertextbook.com7fa cts/2000/ AlinaVayntrub.shtml 
'http://www.stdimension.de/int/Cartography/mwtour.htm 



■^One Hght year represents a distance of approximately 9.5 • lO^^m. 
■^We do not make the distinction between the impact parameter (6) and the minimum 

distance approach (ro). The exact relation 7'o = b ( ^ ^"^^ ^ shows that for small D, one 



has ro ~ b. 
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c, the velocity v{r) is greater than c, the (weak) correction term being of order 
For small D, the minimum distance tq approach is of order b, leading 
to a maximum supraluminal correction for the velocity of order that is 
small. As a temporary conclusion, I would say that for practical purposes, 
the weak gravity limit of the classical case, together with the identification 
^oo c, yields a correct order of magnitude for the (GR) value. Some 
inconsistencies of classical mechanics appear in this identification, but are 
small in the weak gravity limit. 



2.2 The transparent deflector 

A preliminary remark is that D{b) is not a monotonous function, since it 
vanishes both for 6 = and b = oo. For [b > R), the deviation is exactly given 
by equation Q, the weak gravity approximation corresponding to equation 

The case {b < R) requires the study of two distinct phases of the motion, 
since the gravitational field g{r) on the test particle has different expressions, 
depending on the particle position r. We respectively denote these phases 
by (out) and (in). Their respective contribution to the total deviation D{b) 
will be denoted by D(^out) and D(^in). We further assume that the deflecting 
galaxy is homogeneous, and neglect all collisions with the particles in phase 
(in). 

In phase (in) of the motion. Gauss' theorem gives 

r = 9{^n){r) = = -u^r (5) 

This harmonic motion has period Tq = ^ = 271^ The trajectory inside 
the deflector is elliptic or partially elliptic. One should express the initial 
and final boundary conditions to get the deflection D(j„)(6). These boundary 
conditions depend on phase (out) of the trajectory {r > R), where Newton's 
equation reads 

r = 9{out){r) = (6) 
The total deviation D{b) is given by 

D(&) = D(,„)(6) + D(,„i)(6) (7) 
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2.2.1 Phase (out): (r > R) 



(a) Rigorous solution 

We want to calculate the deviation between initial (source S) and final 
(entry into defiector (A)) points. The final point M has position rj„ = OM 
and velocity Vin (Figure 2(a)), with |rj„| = R, and Vin = fo^m- Energy 
conservation yields = v^^ + 2^. The deviation D'^]^ for this part of the 

(out) trajectory is given by cosD^y^f = Ui ■ Um- 

Equating the projections of the Runge-Lenz vector (eq. 0) along Ui for 
the initial (S) and final (M ) points leads to 

QMm = Ui ■ {vin X L - GMm^) (8) 

R 

Defining cos(7r — $) = ■ ^ (Figure 2(a)), we get 

sinD«= 7^(1 -cos $) (9) 

bVoVoo 

The exact solution of the Kepler problem (see eq (14,7), p. 46 of [7]) yields 

- - d -d 
$ = — arccos( f ) + arccos( . ) (10) 
V 1 + (i^ V 1 + (i^ 

with = if . 

(2) 

Taking into account the symmetrical (AE) contribution (deviation DlJ^) 
finally gives 

Doutib) = 2arcsin(-^(l - cos$)) (11) 

bVoVoo 

where $ is given by ()10j). 

(b) Weak gravity approximation 

As previously stated, we expect gravitational deviations to be weak. 
Equation ()11|) shows that the weak gravity ( "small Q" ) limit can be obtained 
by setting d = in eq. (|1U|). The weak gravity approximation therefore reads 



Z,„,(6,.2— (l-Jl--) (12) 



In particular, the small b limit is given by Dout(b) ~ 4, with / 



/' """" J gu 
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2.2.2 Phase (in): (r < R) 



(a) Rigorous solution 



The geometry of the (in) phase is shown in Figure 2 (b). The angle a is 
the angle between the radius vector fin and the velocity Vin as the particle 
enters the deflector (Figure 2(a)). From the conservation of the angular 
momentum, one has 

mvoob = mvoRsiiaa (13) 

where = v^^ + 2^. 

To calculate the deviation Din{b) between the entry and exit points, one 
may solve the harmonic motion of eq. (0). A more convenient way is to use 
the Runge-Lenz matrix invariant associated with the harmonic oscillator 

A = y(cUo'rr + vv) (14) 

Starting from the invariance of A, simple calculations |2j show that 

%„)(&) = 2(a-/?) (15) 

where sina = and tan2/3 = 

(b) Weak gravity approximation 

In this approximation, we obtain sin a ~ ^, and 

Dinib) ~ tan An(&) = tan2(a - /9) ~ -|^sin2a (16) 



oo 



2.2.3 Conclusion on the transparent deflector 

The total deviation D{b), for b < R and in the weak gravity approximation, 
is given by eq. (fT^ and (fTB|). Setting w = -|, we have 

CM 1 - (1 - M^)! 
D{b) = D^,^){b) + D(,„,)(6) = 2 |— ( 4 ^) (17) 

When Voo — > c, this result can be compared to the (GR) result jH IH]. 
Notwithstanding the (ubiquitous) factor of 2 between classical mechanics 
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and (GR), eq. (fT7|) is in agreement with eq. (7) of reference [Hj- Note that 
equation (jl7j) can be also be written in a way similar to eq. (@]), namely 



,2 3 

where M{b) = M(l — (1 — ;^)^) is the partial deflector mass contained in a 
cylinder of radius b. 

For small b, one gets D{b) ~ 3j. The length / = '^^f can be viewed 
as a focal length; its order of magnitude for our model galaxy is 2.5 10^ ly, 
much bigger than Rq ~ 2 10^ ly. Note also that D{b) has a maximum for 
b ~ 0.93i?G, with ~ 1.6 IQ-^ rd. 

Using the results of equations Pjl. (|T7jl . we show in Figure 3 the deviation 
D{b) as a function of the impact parameter b, for opaque and transparent 
deflectors. 



3 Finite distance geometry and multiple tra- 
jectories 

We now use our results to discuss an experimentally more relevant situation, 
where both the source (S) and the observer (E) are at a finite distance from 
the deflector (A). The - Euclidean - geometry is shown in Figure 4. We have 
rs = SqA = SH, te = AE. Given the previous orders of magnitude, the 
angles such as /3 = SEA, 9 = AEY and the deviation D, are assumed to 
be small. This implies in particular that b = AY = r^jtan^ << rE,rs- We 
have from (Euclidean) geometry 

sinSEY _ sinESY 

SY ~ EY ^ ' 

From triangle SHY, one has 

SY^ = SH^ + HY^ = r| + {ve tan 9 - {ve + rs) tan l3)^ (20) 

yielding for small angles, SY ~ r5(l + 0(^^, /9^, 9P)). 
From triangle YEA one has 

EY^ = rl + AY^ = rl{l + tan^ 9) (21) 
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yielding for small angles, EY ~ te^I + 0{6'^)). 

Plugging these values in eq. (fT^ . we get to lowest order in the angles 
6,(3,.. 

SEY ESY SEY + ESY D , , 
~ ~ ~ (22) 

rs Te rs + TE Ts + Te 

Since SEY = e-(3-^-(3,we finally obtain 

Dib)^'-^iJ^-f3) (23) 
rs rE 

Equation (|23p . which expresses the condition that a particle emitted from 
(S) reaches the earth, is represented by the dotted lines in Figure 3. For an 
opaque defiector, one may get one or two solutions for b. For a transparent 
defiector, one may get up to three solutions for b. Rather than studying 
the full problem as a function of /3,rs,rE, we illustrate some particular 
situations 



3.1 A generic case 

This case corresponds to a non zero (3 angle (Figure 4). We focus our in- 
terest on points (1) and (2) of Figure 3, which are the intersections of the 
geometrical equation 

Z,(ft)^!l±l£(A-,3) (24) 

rs rE 

with the b > R gravitational defiection D{b) of Section |21 (see Figure 3). 
Since we have D(b) ~ 2 setting 9 = — leads to 

9^ -(39- 91 = (25) 

where 

rE[rs + rE) 

In this case, there are two trajectories in the S/S.E plane that reach (E). In 
photon language, the observer sees two images {SI) and (5*2) of the source 
(S), on opposite sides of the defiector, with 



^1,2 = ^ = — (27) 

For a transparent deflector, one has another image of the source, correspond- 
ing to the point labeled (3) in Figure 3(b). 
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3.2 Einstein rings 



This case corresponds to the ahgnment of (S), (A) and (E) {j3 = 0). Due 
to the symmetry of revolution around the SAE axis, all trajectories on the 
angular cone = 6e = ^ reach the earth. In photon language, this means 
that the observer sees a ring image of the point source S. For a transparent 
deflector, one also has a direct image. For our model galaxy and r^; ~ r5 ~ 
10^ ly, a typical value is 9e ~ 10~^ rd. 

3.3 The case of a moving deflector 

We briefly consider this case (called microlensing) , because of its experimen- 
tal relevance. Since a detailed comparison with the experiments require the 
use of General Relativity Pfl 12 El? we limit our presentation to orders of 
magnitude calculations. If the deflector (A) moves, with a velocity va, in a 
direction perpendicular to the (SE) axis (Figure 4), the above calculations 
suggest the following scenario: for /5 = (ring image), there is a sudden 
increase in the signal received by the observer, since two trajectories only 
survive for 7^ 0. Physically the transition is gradual, and the observer will 
receive a gravitationally enhanced signal when the position of deflector (A) is 
within a distance 6^; ~ teOe from the full alignment position of the previous 
section ^ The corresponding time interval is Ie ~ For distant sources 
(js » te), an experimental situation corresponding to a sun-like deflector, 
with t>A ~ 200 km s~^, ~ 3.6 10"^ ly and = c, yields an enhanced signal 
during an interval t^; ~ 10® — 10'' s, of order one month. 



4 Time lags 

4.1 A simple optical analogy 

We have seen that several trajectories- or light rays- may reach (E) because of 
the gravitational deflection. Can one further extend the mechanical-optical 
analogy by finding the time lag between the reception of these trajectories 
-or light rays- ? 

simple way to calculate the amplification factor is to use the non linear relation 
9{(3) of equation (|?7jl . The flux emitted by an extended source is proportional to f3d(3, and 
the flux received by the observer is proportional to Odd. The total amplification factor is 
given by ^=||^K + ||fb. 
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We first consider the transparent defiector for b << R, where we found 
D{b) ~ 3 J in section |21 This result may be compared with the optical 
deviation of a spherical glass lens, of radius p and of optical index u which 
reads 

Dopt = - = 7^ (28) 

1^ P Jopt 

where b is the impact parameter of the light ray and fopt = 2(u-i) focal 
distance of the lens. The comparison of the gravitational and optical devi- 
ations suggests that the gravitational deviation may be understood through 
a gravitational index Ugrav-, with ngrav 7^ 1- 

4.2 Particule trajectories and geometrical optics 

The preceding remark can be extended and formalized as follows. Energy 
conservation for a central potential f/(r) reads, in usual polar coordinates 

Tin Tft 
E = —v^ + mU{r) = — (r^ + r^O^) + mU{r) (29) 

Defining tan\l'(r) = we have 

E = ^ . . + mU{r) (30) 

which can be rewritten as 

nu{r) r sin\l>(r) = (- — -)2 (31) 
2mE 

where the "index" nij{r), associated to the potential U{r) is given by nu{r) = 

\/l — Equation is analogous to Bouguer's relation for the propa- 

gation of light rays in a spherically symmetric medium of index nu{r) [TUj . 

We will illustrate this analogy with the case b > R, where U{r) = —2M_ 
The associated gravitational index reads 



, 2gM QM , , 

n,_(r) = Jl + — -^1 + — - (32) 



where we have used the weak gravity approximation. Defining an analog 
Lgrav of the optical path, we may express the time lag 6ti2 between the 
reception on earth of trajectories (1) and (2) of section ItTD as 
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L-2 — Li /(2) ngrav{r) dS2 — /(I) rigravir) dsi 
Oti2 — — [06) 

From eqs. and (j^ . one finds that 6ti2 is the sum of a geometrical part 

d,S2—J dsi 

Staeom = and of a gravitational part 

= !?fi / ;?£i) (34) 

vl^ J(2) r 7(1) r 

The full calculation of the integrals in (see eq. 8.30, p. 240 of [2 ), yields 
a result that depends only logarithmically on the geometrical parameters 
(/3, rg, ...). We therefore estimate St„rav ~ up to a numerical factor 
of order one. Setting v^o = c, we find Stgrav ~ 10^^ s for the sun, and 
St grav ~ 10^ s for our model galaxy. Experiments that confirm the double 
reception of the same "signal" , with a gravitational time lag of order several 
months, can be found in references [H 12 El ISl El ■ 

Finally, it of interest to note that a Fermat approach to the weak gravity 
approximation of (GR) yields an equation similar to equation (jH^ . with a 
(GR) index ~ 1 + 2 2^ HH [121 . 



5 Conclusion 

We have studied, at a qualitative level, a classical mechanical introduction 
to gravitational lens phenomena. This approach rests on the fact that the 
gravitational defiection of a massive particle by a defiector is independent 
of the particle mass. It is only an approximation to the theory of General 
Relativity mHig, but I believe that this "F = ma optics " jlSl CH brings 
together in a very pedagogical way problems of different origins. In particu- 
lar, we have derived orders of magnitude for the weak gravity case, that can 
be compared -up to a factor 2- to the correct (GR) results. As a caveat, we 
have nevertheless pointed out that the particle velocity may become (weakly) 
supraluminal, and this (weak) inconsistency with relativity is to be kept in 
mind. 

Beside the study of General Relativity, the interested student can carry 
further the present approach in several ways. I will only quote here the 
modehng of gravitational lenses by optical lenses of the appropriate shape 
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[H], or the link between equations (0) and (jH)) stemming from conformal 
transformations 

It is a pleasure to thank F. Bernardeau for discussions. 
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Figure Captions 



Figure 1: The deflection geometry for b > R, witli (S) and (E) at infinity. 
Tlie vectors Ui and uj are the unit vectors of the initial (emission) and final 
(reception) directions. The deviation D = arccos(uj ■ Uf) is a function of the 
impact parameter b. The minimum distance approach is Tq, and c::^ b for 
small D. 

Figure 2: The deflection geometry for the transparent deflector 

(a) phase (out) for b < R. The vector Uin is the unit vector along the 
velocity as the particule enters the deflector at point M {vin = voUm). The 
corresponding deviation is given by D^^^Jf- = arccos('Uj ■ Uin) 

(b) phase (in) for b < R. The particle enters deflector (A) at point 
M, and exits at point N. The OX and OY axes are the eigenvectors of the 
Runge-Lenz tensor A. Points M and N are symmetric w.r.t. OY. For clarity 
purposes, the direction of Uin has been rotated with respect to Figure 2(a). 

Figure 3: The qualitative variation of D{b) (a) opaque deflector (b) 
transparent deflector (note the maximum for b ~ Rg)- The dotted lines 
represent various cases of equation ((221) • Solutions (1) and (2) correspond to 
trajectories (1) and (2) of Figure 4. 

Figure 4: A typical geometry for finite distance gravitational deflection. 
Trajectories (1) and (2), deflected by (A), reach the earth (E): the source 
(S) has two images (SI) and (S2). For a transparent deflector, one has a 
third image of the source ( see point (3) in Figure 3(b)), corresponding to a 
trajectory entering (A). 
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